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The quenching of oscillations in interacting systems leads to several unwanted situations, which 
necessitate a suitable remedy to overcome the quenching. In this connection, this work addresses 
a mechanism that can resurrect oscillations in a typical situation. Through both numerical and 
analytical studies, we show the candidate which is capable of resurrecting oscillations is nothing but 
the feedback, the one which is profoundly used in dynamical control and in bio-therapies. Even in 
the case of a rather general system, we demonstrate analytically the applicability of the technique 
over one of the oscillation quenched states called amplitude death state. We also discuss some of 
the features of this mechanism such as adaptability of the technique with the feedback of only a few 
of the oscillators. 

PACS numbers: 05.45 Xt, 87.10 -e, 87.19.Ir 


I. INTRODUCTION 

The interaction among oscillators in a system not only 
leads them into a cooperative dynamics but also of¬ 
ten quenches their oscillations. There are two dynami¬ 
cally different oscillation quenching phenomena which are 
termed as amplitude death (AD) and oscillation death 
(OD) [1, 2j. In AD, the amplitude of oscillation quenches 
to zero, whereas OD is caused by quenching in the fre¬ 
quency of oscillation [1]. It is also defined that AD oc¬ 
curs via the stabilization of a homogeneous steady state 
(HSS) while OD occurs via the stabilization of an inho¬ 
mogeneous steady state The mechanisms underly¬ 

ing these two quenching phenomena have been identified 
recently and the results show that the oarametric mis¬ 
match Q, dynamic 1,1 , time delay |7'hr9] and nonlinear 
couplings [l(| underly AD, while OD occurs mainly be¬ 
cause of the symmetry breaking coupling in the system 
[2]. Recently, diverse routes of transition from the AD to 
OD have also been reported muni- 

Experimental and theoretical studies show definitive 
evidence of oscillation quenching in realisti c sy stems 
ranging from biological [13!, H| , chemical Il5l. Il6j , elec¬ 
tronic [17] and laser d systems to climate [Dj systems. 
Such oscillation quenching in many cases leads to unde¬ 
sirable situations. In the interaction between neuronal 
dynamics and brain metabolism, the decrease in cere¬ 
bral metabolic rate, coupled with the stabilizing proper¬ 
ties of ATP-gated potassium channels, leads to a burst 
suppression in the EEG pattern which symbolizes inac¬ 
tivated brain 0]. This type of suppression results in 
hypothermia, coma and Ohtahara syndrome, a type of 
early infantile encephalopathy and is also observed dur¬ 
ing deep levels of anesthesia. The suppression of nor¬ 
mal sinus rhythm of pacemaker cells causes cardiac arrest 

A- 

Owing to the fatal consequences due to oscillation sup¬ 
pression, interesting efforts have been undertaken to re¬ 
trieve/resurrect oscillations of the system [22, 23j. In [22j . 


the oscillation death in diffusively coupled oscillators has 
been found to be eliminated through a spatial disorder 
in the form of parametric mismatch and in Ref. [23] pro¬ 
cessing delay is used to revoke oscillations successfully in 
delay coupled systems. 

Regarding the above mentioned issues, in this article, 
we demonstrate that the problem can also be well re¬ 
solved by providing a suitable feedback in the system. 
The latter can be found to be present in most of the 
natural systems, including neural networks 0, genetic 
networks [25], vision systems 0, etc. The vital role of 
feedback in controlling the dynamics of the given system 
and the control over synchronization [271 RlOj are already 
known, which can be seen in a variety of fields rang¬ 
ing from electronics [3l|, biology to quantum informa¬ 
tion [32, |33|. For example, the feedback control of deep 
brain simulation has been found to be the most effec¬ 
tive treatment for chronic neural diseases like essential 
tremor, dystonia and Parkinson’s disease [2, [35|. Also, 
the feedback generated by the voltage-gated ion channels 
in neural cells is found to be crucial in generating neural 
signals j36| . 

In this article, we show the applicability of the feedback 
technique in resurrecting oscillations in a wide range of 
systems. We show both numerically and analytically that 
the addition of feedback destabilizes the stable attractors 
which results in a wiping out of the oscillation quenching 
and inducing a resurrection of oscillations. In addition, 
by considering a rather general system, we prove analyt¬ 
ically the above destabilizing nature over the AD state. 

Further, it will be more important to develop an adapt¬ 
able mechanism thereby improving the ones available in 
the literature at present. This is because, for example to 
use the available parametric mismatch method, one needs 
to tune the internal parameters of the system, while the 
processing delay also depends on the underlying process 
of the system where that process may be unknown in 
many situations so as to hinder the efficiency. In con¬ 
trast, the feedback method suggested here can be given 


2 


more easily which is already in practice under different 
contexts such as in deep brain simulation (34], [35| . 

In addition to the above adaptable nature of feedback, 
with the aid of numerical and analytical studies, we show 
the important fact that this method does not impose a 
restriction that the output of all the oscillators need to 
be fedback. From the output of only a few of the system 
oscillators, we show that in typical systems the resurrec¬ 
tion of oscillations can be achieved easily. 

The structure of the paper is as follows. In section 
m we present the general form of the system that we 
consider. In section m we illustrate the role of feedback 
on two important oscillation quenching scenarios, namely 
the symmetry breaking coupling and the parametric mis¬ 
match in a system of diffusively coupled Stuart-Landau 
oscillators, through numerical analysis. In section [IVJ we 
present suitable analytical support of the numerical re¬ 
sults based on an appropriate linear stability analysis. In 
section El we illustrate the role of the considered feed¬ 
back in indirectly coupled or dynamically coupled Stuart- 
Landau oscillators. A realistic chemical oscillator model, 
namely the Brusselator model, is considered in section 
ED We have also proved the applicability of the tech¬ 
nique in more general situations over amplitude death 
state in Appendix [A] In addition, in the Appendix [B] 
we illustrate our method with different coupling schemes 
and with different models such as van der Pol oscillator, 
Rosseler system and so on. Appendix [Cl includes the de¬ 
tails in obtaining the boundary curves of the AD region 
which were given in Sec. IIV1 A summary of our results 
and conclusions are presented in section ivni 


II. THE GENERAL MODEL 

Consider a system of coupled dynamical systems, 

N 

W i = fi(wi) + ey^LjjH(wj) + 7}g(u(t)), 
i=i 

* = 1,2, ...N (1) 

where f) : H d —>> H d characterizes the dynamics of the 
isolated i-th system, w* G R d is a d dimensional state 
vector of the system i, L^ is the d x d coupling matrix of 
the network, (e, rf) G R are respectively the uniform cou¬ 
pling and feedback strengths, H : H d H d is a coupling 
function and g(u (t)) : R d — )> R d is the feedback term 
which can be written as g(u (t)) = Qu(t). Here Q is sim¬ 
ply a d x d constant matrix and u(£) G H d characterizes 
the feedback and it depends on the state vectors of the 
system. Such a dependence of u (t) on the state vectors of 
the system may be linear (Example: u(t) = XlfeLi a fc w /c) 
or nonlinear (Example: u(t) = ^2^=1 a /c( w fe T w/ c ) g w/ c ), 
where a&’s represent weight factors which can take values 
from 0 to 1 and q is a suitable number. In our following 
study, we consider the form of u (t) as u (t) = XlfcLi a k w k- 




FIG. 1: (Color online) (a) Suppression in OD regions of the 
system © with respect to the introduced two different forms 
of feedback (i) — 1 and (ii) ak = &k i, k = 1 , 2 , for the choice 

N = 2 and uj — 2.0. (b)For two choices of e, e = 10.0 and e — 
13.0, the change in the HB point by the increment of P (which 
appears in the expression of a^) with N — 100 in ([2]). Inset in 
(b) projects the OD and oscillatory (OS) regions in the (e, rf) 
space for P — N. (c),(d) and (e): Temporal behaviors of the 
system in the OS, OD and revoked OS state respectively for 
the values of (e, 77 ) = (9.0,0), (14.0,0), and (14.0,0.4). The 
color bar in Figs, (c), (d) and (e) gives information on the 
value of yi corresponding to different colors. 


In the Appendix El we have considered rather general 
forms for and L^ and shown analytically that the triv¬ 
ial AD state which appears in the system could be wiped 
out through the strengthening of 77 so as to resurrect os¬ 
cillations. 


III. DIFFUSIVELY COUPLED SYSTEM: 

NUMERICAL ANALYSIS 

A. Symmetry breaking coupling 

To start with, we use the paradigmatic model known as 
the coupled Stuart-Landau oscillators for the purpose of 
illustration for the advocated feedback method. It is well 
known that the dynamical equation defining the Stuart- 
Landau oscillator can be obtained from a general ordi¬ 
nary differential equation near a Hopf bifurcation point 
[3t m. As the Hopf bifurcation arises widely in the 
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literature, the Stuart-Landau oscillator helps to model 
a variety of systems in different areas ranging from bi¬ 
ology [391-41] to lasers [42 2 43 [ and is also used in the 
reaction-diffusion process 3 a 44]. This model is often 
used in neural networks to model spiking neurons [ 22 t 
l4lj |. The first reason for using the model in neural net- 
works is that the periodically spiking neurons have an 
exponentially stable limit cycle attractor and secondly 
the real part of the complex amplitude of the Stuart- 
Landau oscillator can describe the membrane voltage in 
the neurons and the imaginary part can be related to the 
recovery variable embedding the effects of the other vari¬ 
ables of physiological neuron models ,45]. Thus, in the 
literature we can find the use of this model in studying 
the effects of synchronization and desynchronization in 
neural networks [HI, Hl|, 5(1 - 149) and also various collec¬ 
tive dynamical states, including chimeras [HqI- 152). In our 
study, we also include other useful models such as the 
van der Pol oscillator, Rossler system and Brusselator 
model. The corresponding results are briefly indicated in 
the Appendix B and Sec IVI1 

Now, we first consider a system of coupled Stuart- 
Landau oscillators which is characterized by 


n ( \ I CJiyi Tj^Xi i tj / \ 

WW ‘ ) -| „, + «»- rh, •’ 


L« = sr 


—NSa + 1 0 


0 0 
g(u(t)) = Qu (t), 


u W = F 


Q_ N’ 


C^k^k 


( 2 ) 


where r\ ■= x\ + yf, W[ = [xi yi] T , Sij is the Kronecker 
delta (S^ = 0, if i ^ j and 5^ = 1, if i = j ) and I 
represents the identity matrix, in Eq. (J2]) . Here the dif¬ 
fusive coupling acts only on the first half of the evolution 
equations fa- variable alone) which breaks the rotational 
symmetry (22j and consequently induces oscillation death 
in the system. 

To elucidate clearly the role of feedback in (|2j) , we first 
consider the case N = 2, with uq = co >2 = w. In the case 
6 = 0, r] = 0, we note that the individual systems in (|2]) 
show limit cycle oscillations with fa | = 1. The intro¬ 
duction and strengthening of diffusive coupling (e ^ 0) 
stabilizes the symmetric pair of nontrivial equilibrium 
points {xl,y{,xl,y%) = {a, 

'■* ~ cb\ 2 , b\ a 


n ,2 


± 


i , 6 *, —a*, — 6 *), where i » 1 , 2 , 
and 


1-\-ujc 

1+c 2 


C = 


- £ +vV^ 

2u> ’ 

]} which gives 


through a sub-critical Hopf bifurcation 
rise to OD. With the introduction of feedback 77 7 ^ 0, 
these nontrivial equilibrium points soon lose their stabil¬ 
ity via a super-critical Hopf bifurcation. In our study, we 
introduced such a feedback in two ways: (i) = 1 and 

(ii) ak = Ski, k = 1, 2 (again here S denotes the Kronecker 
delta). First by setting dk = 1, for different values of e we 
traced the Hopf bifurcation points and these points are 
collectively shown as the HB curve—1 in Fig. HJa). The 
region lying under this curve is an OD region (denoted 
by OD— 1 ) and the parametric region above this curve is 


free from OD and corresponds to oscillatory states (OS). 
Similarly for = <5^i, k = 1,2, the OD region (OD-2 
which includes OD-1 also) and the curve of Hopf bifur¬ 
cation points (HB curve—2) have been shown in figfflfa), 
which show that the uniform distribution of ak as in (i) 
helps to redeem from the OD state sooner than in the 
case (ii). 

By extending the constituents of the network to N = 
100, we have verified that this technique can work as well 
with larger N. Further, we have checked whether the 
feedback needs contributions from all the constituents 
of the network. This is vital as in a practical situation 
we cannot assure or impose all the constituents to con¬ 
tribute to the feedback. Thus, we have distributed a/ds 
as dk = , & = 1,2, ...,7V, where P determines 

the number of contributing components of the network. 
By varying P, we have drawn the HB curves separating 
the OD state with the OS state for two different values 
of e, e = 10.0 and e = 13.0 in Fig. QJb). Interestingly, 
these curves demonstrate clearly that the contribution 
from even a single oscillator is sufficient to revoke oscil¬ 
lations in the network. Secondly, the critical value of 
r] (HB point) above which the OS state arises gets de¬ 
creased sharply with that of P. These facts prove that 
this technique can work well regardless of the number of 
oscillators present in the network and the number of them 
which contributes towards the feedback. For simplicity, 
we chose P = N in the following studies. 

The inset of Fig. []Jb) depicts the information about 
OD and OS regions with N = 100 oscillators. Corre¬ 
spondingly, in Figs. mMd) andQJe), we have captured 
the temporal behaviors of the system (in the y- variables) 
for different sets of (e, 77 ), for a finite time interval af¬ 
ter leaving out sufficiently large transients. The first one 
(Fig. [He)) shows the temporal behavior at an OS state 
of the system when 77 = 0 and e = 9.0, where the value of 
e is not sufficient to induce OD. Now increasing e to 14.0 
(while keeping 77 = 0), the subsequent figure (Fig. 0Jd)) 
shows the quenching of this oscillation. Now switching 77 
on, the temporal behavior in Fig. []Je) shows the resur¬ 
rected oscillations for 77 = 0.4. 


B. Effect of parametric mismatch 

Rubchinsky and Sushchik |22] have introduced a dis¬ 
order in the form of parametric mismatch which revokes 
oscillations in (| 2 ]) . On the other hand, just like symmetry 
breaking which is predominant in inducing OD, the para¬ 
metric mismatch is also a key candidate that induces AD 
in the system. Recently, Koseska et al. in 0, [53|, [H| have 
shown that an increase in this inhomogeneity not only in¬ 
duces AD but also OD, whereas the feedback mechanism 
that we consider here does not show such a behavior in 
the system. This feature provides a definitive advantage 
over parametric mismatch. 

Now, we augment the system with a parametric mis¬ 
match. For the mismatch in the parameter cj, — = 
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FIG. 2 : (Color online) The transition from AD (a) or OD (b) 
state to OS state via Hopf bifurcation with respect to 77 in the 
system © endowed with a parametric mismatch of A = 5.0, 
UJ 2 — 2.0 and N = 2. (c) The reduction of AD regions in the 
(A, e) space with the increase of 77 and by keeping U 2 — 2 and 
N = 2 . (d) and (e) The emergence of OS state from AD and 
OD states in N = 100 case for A = 5.0 and 0 J 2 = 2.0 which is 
obtained for two different set of initial conditions. AD region 
in Fig. (e) is also in conformity with the analytical results 
presented in the Sec. m 


A = 5, in the above coupled Stuart-Landau oscillators 
0-0 with N = 2, we find the existence of both AD 
(e = 6.0) and OD (e = 10.0) while 77 = 0. Now, from 
these death states the transitions towards OS state by 77 
are demonstrated in Figs. 0a) and 0b), which show the 
destabilization of both the AD and OD states via super¬ 
critical Hopf bifurcations. Further, the role of 77 over A 
and e is more clearly demonstrated in Fig. 0c), where 
the colored islands denote the AD regions for the values 
of 77 = 0 . 0 , 0.3, 1.0 and 1 . 8 . One can check that similar 
phenomenon occurs to OD regions also which we do not 
depict here explicitly. 

Now considering the case of N = 100 globally coupled 
oscillators with uji = 10.0 for i = 1, 2,..., 50 and uji = 2.0 
for i = 51, 52,..., 100 in 0, we have demonstrated the 
reduction in the AD and OD regions in Figs. 0d) and 
He). The above two figures are plotted for two different 
sets of initial conditions. Among them, the AD region 
in Fig. 0e) is the analytically relevant region which has 
been obtained in the next section. 

In the next section, the above obtained numerical re¬ 
sults on the system m are verified through analytical 
results wherever possible. Also, we illustrate the ap¬ 
plicability of the technique to several situations in Ap¬ 


FIG. 3: (Color online) Without parametric mismatch case: 
Fig. (a) Analytical results of the OD regions in (e,rj) space 
for N = 2 oscillators in system © with cji = 0 J 2 = 2.0. With 
parametric mismatch case: Fig. (b) AD regions of system 0 
for A = 4, UJ 2 — 2 and a\ — (12 = 1 bounded by the curves 
defined in (1151) . Fig. (c) Analytically obtained boundaries of 
AD regions in (A, e) for different values of 77 (=0, 0.3,1.0,1.8), 
where a\ — <22 = 1. Fig. (d) HB curves for the three different 
cases (i) ai = 02 = 1 , (ii) ai = 1 , 02 = 0 and (iii) a\ — 
0, <22 = 1 for A = 4 and UJ 2 — 2.0. 


pendix[Bj where we considered (i) Repulsive link: Stuart- 
Landau (SL) oscillator, (ii) Conjugate coupling: SL os¬ 
cillator, (iii) Repulsive link: van der Pol oscillator, (iv) 
Directly and indirectly coupled Rossler system and (v) 
Other chaotic oscillators (Sprott and Lorenz systems). 


IV. ANALYTICAL CONFIRMATION OF 
SUPPRESSION OF DEATH STATES 

In this section, we present relevant analytical confirma¬ 
tions of the numerical results corresponding to the system 
0 for the cases with and without parametric mismatch 
discussed earlier in Sec. m In addition through the ob¬ 
tained analytical results, we show the effectiveness of the 
technique with the feedback contribution coming from a 
few number of oscillators in the network. 


A. Without parametric mismatch: N = 2 case 

In the absence of any parametric mismatch, as pointed 
out earlier the system in 0 has a trivial equilibrium 
point Eq: (0, 0, 0, 0) and two pairs of non-trivial equilib- 
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rium points as 

Ei : (xl,yi,xl,yZ) = (a*,b*,-a*,-b*), i = 1 , 2 ,3,4, 


where 



a l,2 = c ^l,2 5 ^1,2 = ±\J 

ll+UJC 


1 + c 2 ’ 


<4 = db* 3A ,bl 4 = ±y 

1 IT ujd 

(3) 

i + d? ’ 

where e = 

and d = 

The 


eigenvalues of the Jacobian matrix 


J = 


1 _L H 

2 ' 2 


/ Ai —uj — 2 x*yi 
c j — 2x*y* B\ 0 

0 A 2 —id — 2 x^V2 

0 f u-2x* 2 yl B 2 


£ + 2? 
2 ' 2 


(4) 


where 1 — 3x * 2 —?/* 2 —+ B\ 1 —x * 2 —+ 

i = 1 , 2 . While 77 = 0, the eigenvalues of J corresponding 
to the trivial equilibrium point Eq are 


d °2 = h 2 - e ± \/e 2 - 4w 2 ), /4°1 = 


= 1 =b icj. 


(5) 


The eigenvalues corresponding to the equilibrium points 
Ei and E 2 are given by 




= -sJ7> -4u> 2 , = (e - 2) - Ve 2 -4w 2 , 

/4m = ((e - 1) - \/e 2 - 4w 2 ) ± VCl , j = 1, 2 (6) 


where 


Ci = 


—2(e — 1) + \/C"^4cj2(y / e2"TT4j2 _ ( e _ 2)) 


(7) 


Similarly, £3 and £4 have the set of eigenvalues 
/i^ — \/ e 2 — 4a; 2 , /i^ = (e — 2) + \/e 2 — 4oj 2 , 


= ((e - 1) + G 2 - 4^ 2 ) ± VC~2- j = 3,4, (8) 

where 

—2(e - 1) + Ve 2 - 4a; 2 (Ve 2 - 4cc; 2 + (e - 2)) 

0 2 = - - -(9J 

From Eqs. ( 15151 ) we can note that among the five 
equilibrium points, E\ and E 2 are found to have all 
their eigenvalues satisfying the condition Re[p\ < 0 

for the parametric range e > ( 1+4uj ) while uj > 

1+4uj2 < e < | (2 — |Vl — 3 cj 2 ) while | < uj < 


1 

3 ’ 
| and 


2 3 V 3 v 2 

e > | (2 + | y/1 — 3cj 2 ) while cj < and thus they are 
stable in this range. On the other hand the other equilib¬ 
rium points Co, £3 and £4 can never become stable for 
any choice of parametric values. Thus the stabilization 
of Ei and E 2 essentially gives rise to oscillation death in 
the system. 

Now we introduce the feedback in such a way that 
ai and d 2 of u (t) in (| 2 ]) take the values ai = 1 and 
a 2 = 1. With such a choice, we find that the stability 


determining eigenvalues corresponding to the equilibrium 
points change as 

= and ^4 = Mm + V, j = 0,1, 2,3,4. (10) 

The above equation shows that an increase in 77 can desta¬ 
bilize the equilibrium points Ei and E 2 through a Hopf 
bifurcation. In the cases where uj > 1, for all values of e 
the Hopf bifurcation occurs at 

77 = (1 - e) + \/e 2 - 4cj 2 , (11) 

whereas in the case of uj < 1, if e > ^ 2uJ , the Hopf 

bifurcation occurs at 

v = (1 - e) + \/e 2 - 4 w 2 - y/c[ ( 12 ) 

where C\ is given in ( 171 ). If e < — 1 1 2 ( ))' +1 ' ) , the Hopf bi- 
furcation occurs as in CD- The other equilibrium points 
£o, £3 and E 4 are found to remain unstable. For the 
case of uj = 2.0, the curve of Hopf bifurcation points 
(77 = 1 — e + \/e 2 — 4cj 2 ) separating the death regimes 
with the oscillatory regimes is shown in Fig. [3](a) which 
matches exactly with the one obtained numerically (HB 
curve - 1 in Fig. DJa)). The analytical treatment of the 
other case corresponding to the HB curve-2 can also be 
investigated in a similar manner, though the results can¬ 
not be written down in such a transparent manner. So 
we do not present the details here. 


B. With parametric mismatch: N = 2 case 

Next with the introduction of a parametric mismatch 
in the system ©, we prove the validity of the technique 
for the trivial AD state analytically for the N = 2 case. 
As mentioned earlier in Sec. mm an increase in the 
value of the parametric mismatch parameter A causes 
the stabilization of the trivial equilibrium point (0, 0, 0, 0) 
and thus introduces AD in the system. To destabilize the 
latter, we first introduce the feedback in such a way that 
g(u(t)) - (£), where u(t) = ai wq + a 2 w 2 . 

As before, through a linear stability analysis of the 
system ([ 2 ]) with parametric mismatch included, we look 
for the stable regions of the equilibrium point (0, 0, 0, 0), 
which can be studied from the characteristic equation of 
the linear eigenvalue problem of the system as 

/i 4 + A 3 /i 3 + A 2 /i 2 +■ Aiii + Aq = 0. (13) 

The coefficients A 3 , A 2 , Ai and Ao in the above equation 
are given by 

A 3 — e — 4 — J 77 , 

A 2 = ^a 2 77 2 + <277(3 — e) + (6 — 3e + + cj 2 ), 

Ai = A 1177 2 + Ai 2 77 + A 13 , 

Ao = Aoirj 2 + A0277 + A03, 
where 
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time t 


Ci 2 = 1 (only the low frequency oscillator contributing). 
The analytically obtained Hopf bifurcation curves for all 
the above three cases have been presented in Fig. [3jd). 
From the figure, we can note that for the first two cases (i) 
and (ii) the OS state gets revoked from the AD state for 
even small values of 77 and the Hopf bifurcation curves of 
these two cases are closer to each other. But, in the case 
where the low frequency oscillator alone is contributing, 
we find that comparatively higher values of 77 are needed 
to revoke oscillations. This shows that for a quicker res¬ 
urrection of oscillations, the feedback from the high fre¬ 
quency oscillator is preferable. However, the resurrection 
of oscillations is possible even if one of the oscillators is 
contributing towards the feedback. 


FIG. 4: (Color online) Temporal behaviors of the system J2} 
with N = 100 oscillators in the OS, OD and revoked OS states 
respectively, for the values of (e, 77 ) = (4.0, 0.0), (7.0, 0.0), and 
(7.0, 3.0). Here, we have taken N± = 50, A 2 = 50, pi = 20, 
P 2 = 25, LJ 2 = 2.0 and ui = 10.0. The color bars in Figs, (a), 
(b) and (c) represent the values of the variables yi. 


An i(e —2 )a 2 , 

A 12 — (2e — 3)n — (c<qa 2 0J2 a 1)5 
A13 = (3e-4) + i(e-4)(cj? + <*;£), 

A 0 i = i [(1 - e)a 2 + (cjia 2 + a; 2 ai) 2 ] , 

A02 — i [(4(1 - e) + aJiUJ2)d - (e - 4)(a;? a 2 + co^ai)] , 
A03 — (1 + CJi )(1 + CJ 2 ) — — (2 + + CJ 2 ). (14) 

Here a = a\ + < 22 . Since the characteristic equation 
for the eigenvalues is quartic in nature, we use the well 
known Routh-Hurwitz criteria [55j to obtain the stable 
AD regions of the system. By doing so, we find that the 
AD regions are bounded by the curves 


and 


Aq = 0 or r\ — 


-A 02 ± ^/Ap 


02 


4AoiAq3 


2 A, 


01 


A 3 A 2 A 1 — Af — A 3 A 0 — 0. (15) 


The details of obtaining the boundary curves from the R- 
H criteria are presented in the Appendix [CJ and the AD 
region bounded by the curves (H5|) has been shown in Fig. 
[3l(b) for a\ = 1 and <22 = 1. Fig. [3](c) which portrays the 
boundaries of the AD regions in the (A,e) space clearly 
shows that the analytical results match nicely with that 
of the numerical results given in Fig. He). 

Next by varying the nature of oscillators contribut¬ 
ing towards feedback we have plotted Fig. Etd), where 
we considered three cases (i) a\ = 1 , <22 = 1 (both the 
oscillators contributing) (ii) a\ = 1 , <22 = 0 (only the 
high frequency oscillator contributing) and (iii) a\ = 0 , 


C. With parametric mismatch: N = 100 case 


Next, we extend our studies on the revocation of oscil¬ 
lations from the AD state for the case of N = 100 oscil¬ 
lators, where the parametric mismatch in the system is 
introduced in such a way that the system has two groups 
of oscillators. The first group contains N± oscillators with 
uji = uj 1 , i = 1,2,..., TVi and the second group has N 2 os¬ 
cillators with uji = CJ 2 , i = 1, 2,..., AT 2 with > UJ 2 and 
Ah + Ah = N. Also, we consider that among the Ah os¬ 
cillators in the first group only the output of a sub-group 
of pi oscillators is fed back, in other words, the a/ds of 
u (t) in (|2j) take the values as 

a k = 1, for k = 1,2, ...,pi, 

a k = 0, for fc =pi + l,pi + 2,Ah. (16) 

Similarly in the second group of Ah oscillators only the 
output of P 2 oscillators is fedback or 


a k — 1, for k — Ah + 1 , Ah + 2 ,..., Ah + P 2 ? 
a k = 0 , for fc = Ah + P 2 + 1? Ah + £>2 + 2, •••, N{17) 


The total number of oscillators contributing towards 
feedback is N p = p\ + P 2 . 

The results corresponding to N\ = N 2 = 50 oscillators 
presented earlier in Sec. IIII Bl clearly demonstrate the 
appearance of AD state in the system. The temporal 
behavior of the system in the original OS state, AD state 
and revoked OS state are shown in Fig. SJ which shows 
the coherent nature among the oscillators in the first and 
second group. Next to study the case of the N = 100 
oscillators analytically, we first try to reduce the problem 
to a simpler level. 

Due to the existence of coherence among the oscillators 
in the first and second group, we represent the state of 
the oscillators in the first group by z\ and the state of 
the oscillators in the second group by Z 2 


fi(Zi) 


Xi - ujiYi - RfXi \ 
Yi + LJiXi - R\Yi ) ’ 


h ( z j) = Z P 


/(_2^. + l)(l-p) 0 \ 

v o o;’ 


2 

u(t) = y^afcZfc 

k =1 


g(u (t)) = Qu (t), 


Q = I. 


(18) 
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FIG. 5: (Color online) Fig. (a) shows the numerical (contin¬ 
uous line) and analytical (dotted line) results on the bound¬ 
aries of AD regions in (e, 77 ) space for the values of (pi,p 2 ) = 
(50, 50), (20, 25) and (10,15). Fig. (b) shows the curves of rj c 
for various combinations of (pi, p 2 ) with fixed N p = pi +£>2 = 
45 and e — 7.0. Fig. (c) Fixing pi — 50 or p 2 = 50, the 
values of r) c for different values of p 2 or pi is plotted with 
e — 7.0. Fig. (d) shows the value of r\ c for different values of 
pi and p 2 for e = 7, where the values of r) c are represented by 
the color function. In the above figures, we make the choice 
Ni — N 2 = 50, uji m 10.0 and UJ 2 — 2.0 


where Zi = [Xi Yi] T and P 2 = X? + if, i = 1,2. The 
coefficients in the feedback take the form a\ = ap and 
a 2 = ( 3(1 — p), where p = , which is the ratio of oscil¬ 

lators present in the first group, and 1 — p = ^ is the 
ratio of oscillators present in the second group, a and (3 
are the ratio of the oscillators contributing towards feed¬ 
back from the first and second groups, respectively. Thus 
a = and (3 = J^-, where p\ and P 2 are respectively 
the number of oscillators in the first and second groups 
that are contributing towards feedback. The stability of 
the system (fl8l) around the equilibrium point (0,0, 0,0) 
can be studied similar to the previous case and one can 
obtain the characteristic equation for the eigenvalues of 
the above system m as 

lP + B 3 fp + B 2 P 2 + B\i± + Bq = 0 . ( 19 ) 

The coefficients in the above equation are given by 

B 3 = (e — 4) — 2Pr), 

B 2 = P 2 r ] 2 - 2(e - 3)Pr] + (6 - 3e) + 

Bi = Burj 2 + B 12 rj + B 13 , 

Bo = Boip 2 + B02V + Bo 3 , 
where 

Bll = (e-2 )P 2 , 


B 12 — 2/3(1 — p)(— 3 + 2e — oj±) + 2ap(—3 + 2e — cj|), 

B 13 = (3e + 2) — 2 (uj\ + cj|) + e(puj \ + (1 — p)cu 2 ), 

B 01 = (1 - e)P 2 + (pau 2 + (1 - p)/3u 1 ) 2 , 

Bo 2 = 2P(1 - e) + 2(/3(l — p)u\ + paul) 

+ep( 1 - p)(uj 1 - uo 2 )((3uoi - <acj 2 ), 

B 03 = (1 + )(1 + cu$) — e (l + puj\ + (1 — p)idl) (20) 

where P = (pa + (1 — p)/3). As in the previous case, 
using the R-H criteria we determine the AD regions of 
the system. Consequently, we find that the AD regions 
are bounded by the curves defined by 

n n —B 02 ± \/ Bq 2 — 4B 01 Bo3 

Bo = 0 or r] = -- 

znoi 

and 

B 3 B 2 B x -B\- BjBo = 0 (21) 

Next, using the above relations we find the AD regions in 
the different cases of the system. We consider here only 
two cases, (i) N\ = A 2 , (ii) N\ / N 2 . 

Case-1: N\ = N 2 = 50. In this case, the population 
of the high frequency oscillators (N\) and the popula¬ 
tion of the low frequency oscillators (7V 2 ) are equal (Note 
that u)\ > cj 2 ). First, we check the consistency of the ob¬ 
tained analytical results with the numerical results. For 
N\ = N 2 = 50, we have plotted the boundaries of the 
AD regions obtained from analytical and numerical stud¬ 
ies for different values of p\ and p 2 in Fig. [5fa). (Note 
that for pi = N\ = 50 and p 2 = /V 2 = 50, the numerical 
results have been given in Fig. E(e) of Sec. nng. The 
latter shows the consistency between the numerical and 
analytical results. 

Now, we look at the preferential feedback configuration 
for quicker resurrection of oscillations. For the purpose, 
we fix the total number of oscillators contributing to¬ 
wards feedback as N p = 45 and vary the number of oscil¬ 
lators contributing from the first group (p \) and from the 
second group (p 2 ). Fig. [Sfb) shows the boundaries of the 
AD regions for different values of p\ (also p 2 = N p — pi), 
where we can observe that on increasing the contributions 
from the high frequency oscillators (or pi), the resurrec¬ 
tion of oscillations occur for lower values of 77 . Again in 
Fig. Etc), we first fixed p\ = 50 and plotted the critical 
value of r] needed for resurrection of oscillations (r] c ) for 
different values of p 2 . Similarly, we fixed p 2 = 50 and 
plotted r\ c for different values of p\ in the same Fig. [5jc). 
From the figure, we find that a considerable decrease in 
the value of r\ c occurs only when pi is varied. The above 
results show that a feedback from the high frequency os¬ 
cillators is more preferable for a quicker resurrection of 
oscillations. They are also evident from the Fig. [5fd), 
where r) c is represented by a color function. In the figure, 
one can find that the rate of change in the value of rj c is 
larger along the p\ direction. 

Case-2: N± 7 ^ /V 2 . Just as the technique shows pref¬ 
erence over the feedback of high frequency oscillators, it 
shows dependence over the population in the two groups 
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60 


(b) 


40 

20 


10 20 30 40 


P i 


Pi 


FIG. 6 : (Color online) The critical values of r\ c for different 
values of (pi,p 2 ) for the case (a) Ni > N 2 (N± = 60, N 2 — 40), 
(b) N± < N 2 (Ni — 40, N 2 — 60), with uj\ — 10.0 and 
u 2 = 2 . 0 . 


of oscillators. To illustrate the above, we considered two 
situations (i) Ni > N 2 , (ii) Ni < N 2 . Fig. ( 6 ^a) has been 
plotted for the case (i) where Ni = 60 and N 2 = 40. 
Similarly Fig. ETb) corresponds to the case (ii) where 
Ni = 40 and N 2 = 60. The values of r] c for different 
values of pi and p 2 are presented in Fig. [ 6 j Compar¬ 
ing Fig. Ha) with Fig. [Hb), we find that the values of 
p c are larger in the case of N± < N 2 than that in the 
case Ni > N 2 . Thus in the case where lower frequency 
oscillators are highly populated than the high frequency 
oscillators, we require stronger feedback to revoke oscil¬ 
lations. 


7 c 

14 

10 

6 

2 





t 100 
o 

o 

o 50 


e = 12.0, 77 -O.O 


5 

time t 



FIG. 7: (Color online) (a) The OD and AD regions in the 
(e, 77 ) space corresponding to the system (1221) with N = 2 and 
uj — 3.0, (b) The reduction in the AD and OD regions in the 
system with N — 100 coupled oscillators in the (e, 77 ) space 
for uj — 3.0. (c), (d), (e) and (f) The corresponding temporal 
behavior of the system (I22D of N — 100 oscillators in the AD, 
OD and the revoked oscillatory state. 


V. DYNAMICALLY COUPLED SYSTEMS 

Following the studies on directly coupled systems, we 
turn to check the validity of the proposed scheme to indi¬ 
rectly coupled oscillatory systems. For this purpose, we 
consider a collective system coupled to a dynamic envi¬ 
ronment, where the combined set of dynamical equations 
is characterized by 


( Xi - ujyi - r^Xi \ 

Hi + UXi ~ TiPi J , 

J 


H(wj) = wj, L ij = 


g(u (£)) = Qu(t), Q = — 




( 22 ) 


Here, the state of the system along with the environ¬ 
ment is defined by the state vector w 7 = [xi yi ^] T , 
where the variables Xi and yi correspond to the system 
and vi represent the environment. While 77 = 0, the 
increase in e causes a stabilization of the trivial equi¬ 
librium point (x*,y*,v*)=( 0 , 0 , 0 ) and a further increase 


in e stabilizes a pair of nontrivial equilibrium points de¬ 
fined by (x*, 7 /*,<)=((-l)^*,(-l)^*,(-l)V), where 

a*= kb*, b* = ±^$0, c* = and k = £(-3e + 

V9e 2 — 16cj 2 ) via pitchfork bifurcation 0 . 

The introduction of the feedback, 77 7 ^ 0, destabilizes 
both the AD and OD states via Hopf bifurcation. Such 
a reduction in the territories of AD and OD states is 
illustrated in the (e, 77 ) space in FigJTJa). The curve made 
up of the Hopf bifurcation points separates out the AD 
and OD regions with the OS region. 

By extending N to 100, Fig. [71(b) depicts the AD and 
OD regions of the system. When 77 = 0, the AD state 
which arises by increasing e is found to disappear with an 
increase of e. But for larger e, we find the appearance of 
the OD state. Now by switching 77 on, both the AD and 
OD states are shown to be wiped out simultaneously. 
The temporal behaviors of the system in the AD state 
and the resurrected oscillatory state which arise through 
the enhancement of 77 are shown in Figs. [71(c) and[7Jd). 
Similarly, Figs. [7Je) and m show the behavior of the 
system at the OD state and the resurrection of oscilla¬ 
tions by an increase in 77 . 
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FIG. 8 : (Color online) (a) The existence of stable limit cycle 
oscillations in the isolated Brusselator oscillator (l23l) for A = 2 
and B = 10. (b) The occurrence of OD through the coupling 
of the system for e = 0.5 and 77 = 0, and (c) The resurrection 
of oscillations by introducing feedback 77 = 1 . 


VI. COUPLED BRUSSELATOR OSCILLATORS 

In this section, we consider an interesting coupled 
chemical oscillator modeled by the Brusselator model 
[56j . We consider the case where two identical cells are 
coupled. In such a case, the functions in Eq. describ¬ 
ing the system are given by 


whether it can control oscillation quenching tendencies. 
For this purpose, we have demonstrated the effect of feed¬ 
back over quenching induced by parametric mismatch 
and symmetry breaking (the key candidates for induc¬ 
ing AD and OD), indirect coupling and for some more 
cases (given in Appendix [Bj , through numerical as well 
through analytical studies wherever possible. 

Further, through analytical studies on AD state of a 
more general system (Appendix [A]), we found the general 
applicability of the mechanism where the feedback resur¬ 
rects oscillations from the AD state. In the case of OD, 
a proper/suitable form of linear feedback would be help¬ 
ful to resurrect oscillations or one can also explore the 
role of nonlinear feedback in such nontrivial OD states. 
From the results obtained for different cases, we find that 
the trivial AD state is found to be destabilized through 
Hopf bifurcation whereas the nontrivial OD states are 
found to be destabilized even through saddle node type 
bifurcation (see Appendix [0. 

In addition to the adaptability of the technique in prac¬ 
tical situations, we have illustrated here one more impor¬ 
tant feature of the technique namely that it does not put 
any restriction over the number of oscillators contribut¬ 
ing towards feedback. Even with the feedback from a 
few number of oscillators we can break the death state 
of the system and thus provide an attractive methodol¬ 
ogy in practical situations. Considering a two population 
network, the contribution from the high frequency oscil¬ 
lators are found to be more preferable compared to the 
feedback from the low frequency oscillators. 
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f»(w») 


L ij — 


— {B + 1 )xi + Xiyi + A 

Bxi - x\yi 


H(wj) = wj, 


-2 Sa + 1 


0 


0 —2 5ij + 1 

g(u (*)) = Qu (t), 


i(*) = E 


afcWfc 


Q = 


In the absence of the coupling, the system shows stable 
limit cycle oscillations which has been illustrated in Fig. 
Ufa). By introducing the coupling, the system tends to 
an inhomogeneous steady state [56j. For example, for 
A = 2, B = 10 and e = 0.5, the system tends to an 
inhomogeneous steady state which has been illustrated 
in Fig. [Sfb). In such a realistic example, by introducing 
the feedback, we found that the oscillations are revoked 
by increasing 77 and it has been illustrated in Fig. |SJc). 
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Appendix A: Destabilization of AD in a general 
model 


In this appendix, we show the applicability of the feed¬ 
back technique over the AD state of a general two coupled 
system. For this purpose, we assume N = 2 and the gen¬ 
eral forms for fp H and Lp in Eq. ®, where can be 
chosen as a polynomial in wp 


VII. CONCLUSION 

From a knowledge of the role of feedback in controlling 
the dynamics and the coherent activities of the system 
such as synchronization p7l - [29j |. we have here analyzed 


fi(Wi) 


Fi(xi,yi) +aiXi + biyi \ _ / Xi 

Gi{xi , y%) T CiXi T diyi J y y% 


i — 1 , 2 . (Al) 


In the above, F{ and Gi are nonlinear functions in aq 
and yi and the constants eq, 6 p eq and di are system 
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parameters. The function H(wy) can be written as 

H(wj) = wj + H(wj), (A2) 

where H is a nonlinear function in w j that introduces 
nonlinear coupling in the system. The coupling matrix 
L ij can be taken as 

where k\ = (—2 Sij +1), and k 2 = + !)• 

The systems are coupled through both direct and con¬ 
jugate variables, where k\ introduces direct coupling, k 2 
and k 2 introduce conjugate coupling, e, a and (3 are 
coupling strengths. The system has a trivial equilibrium 
point at (0,0, 0, 0), which may become stable due to the 
parametric mismatch in the system or due to the cou¬ 
pling in the system. For example, in mi the conjugate 
coupling in the system induces AD even when the os¬ 
cillators are identical. The feedback can be given as 

u(t) = J2k=i w /c, where one can also add nonlinear terms 
in the feedback, if needed: 

g(u(*)) = (A4) 

First, considering the case of coupled identical oscilla¬ 
tors di = a, bi = b, Ci = c and di = d (i = 1, 2), the lin¬ 
earization around the trivial equilibrium point (0, 0, 0, 0) 
can be done. The eigenvalues corresponding to the case 
can be obtained easily (note that the nonlinear terms 
in CD and (jA2|) do not play any role in the linearized 
equation for the trivial equilibrium point). The obtained 
eigenvalues corresponding to the case are of the form 

Mi,2 = Mi,2 + M (A5) 

M3,4 = M3,4- (A6) 

In the above pq :2 and Jl 3 ^ are eigenvalues corresponding 
to the trivial equilibrium point of the system when the 
feedback is absent (77 = 0). They are given by 

Mi ,2 = ^(a + d- (1 + P)e) 

V( a ~ d + (/3 — l)e) 2 + 4(/3e + c)(b + e),(A7) 

M 3,4 = — (a + d — 2(1 + cd)e — (1 + /3)e) 

zb^/(a- — d + 2(o — l)e (/3 — 1)e) 2 H- 4C7, (A8) 

where 

C — (c — /3e)(b — e ) (A9) 

Depending on the values of the system and coupling 
parameters, the real part of the eigenvalues /ip 2 and /i 3,4 
in {SB and (1A6|) are positive or negative while 77 = 0 . 
When all the eigenvalues in (jA5[) and (1A6[) have negative 
real parts, the stabilization of the equilibrium point gives 
rise to AD in the system. From Eq. (IMl) . we notice that 


the increase in 77 causes the eigenvalues /ip 2 to be more 
positive. Thus a destabilization of the equilibrium point 
(0,0, 0,0) occurs or it wipes off AD. If the equilibrium 
point was unstable while 77 = 0 , the increase in 77 never 
stabilizes the equilibrium point. Thus the above analysis 
makes clear the role of 77 in destabilizing the attractor at 
( 0 , 0 , 0 , 0 ). 

We can observe a similar effect even in the case where 
parametric mismatch is present (the case where a\ 7 ^ 2 , 
/ & 2 , ci / C 2 and d± 7 ^ 6 / 2 ) in the system. The 
eigenvalues of the system can be obtained by solving the 
equation 

M 4 T A. 3 /i 3 + A 2 M 2 T A\fi + Aq = 0, (A10) 

where A 3 = —(ai +a 2 ) - (<A +d 2 ) + 2 e(l + a) + 2 e(l + /3) - 
277 . As A 2 , A 1 , Ao are not simple in their form, we do not 
present them here. Although the eigenvalues fi obtained 
from (|A 10 |) are not of the simple form, the stability of the 
equilibrium point ( 0 , 0 , 0 , 0 ) in the different parametric 
regions can be found through the Routh-Hurwitz (R-H) 
criteria. From the R-H criteria, an equilibrium point is 
said to be stable only when all the conditions given below 
are satisfied by the coefficients in the eigenvalue equation 
(ETUI) . The R-H criteria are given as 

Ai > 0 , i = 0 , 1 , 2 , 3 , 

A 3 A 2 - Ai > 0, 

A 3 A 2 A 1 - Aj- AjA 0 > 0. (All) 

If the coefficients in the characteristic eigenvalue equa¬ 
tion (1A10[) fail to satisfy any one of the condition given 
above, the equilibrium point becomes unstable. In this 
aspect, we consider one of the simpler condition in (jAll[) , 
namely A 3 > 0. The condition A 3 > 0 is broken when 
77 > \ (—(cl\ + a 2 ) — (d\ + d 2 ) + 2 e(l + cr) + 2 e(l + /?)), 
thus this clearly shows that an increase in 77 destabilizes 
the equilibrium point (0, 0, 0, 0). Further, more clear ana¬ 
lytical illustration on the role of 77 in the parameter mis¬ 
matched system is given in Sec. IIVBI and IIV Cl with 
Stuart-Landau model as an example. 

As the above type of proof for the non-trivial OD state 
is too cumbersome, we have illustrated the role of feed¬ 
back over the state with more examples in the body of 
the paper as well in Appendix [B] both numerically and 
analytically (in some cases). From the above illustra¬ 
tions, one can also notice that in the case of AD, the 
nonlinear feedback terms cannot play any role (as they 
lose their significance in the linearized limit) and do not 
provide any control over it, whereas in the case of OD, 
the nonlinear feedback also can provide a control over it. 

Appendix B: Additional Examples 
1. Repulsive link 

We consider the case of two Stuart-Landau oscillators 
coupled diffusively with a repulsive link (N = 2), as stud¬ 
ied in [57]. The functions characterizing this equation 
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have the forms, 
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FIG. 10: (Color online) Temporal behavior of the system (IBID 
in the (a) AD state for (e, 77 ) = (3.0, 0.0), (b) the revoked OS 
state for (e, 77 ) = (3.0, 3.0), (c) OD state for (e, 77 ) = (4.5, 0.0), 
and (d) the revoked OS state for (e, 77 ) = (4.5,12.5). 


FIG. 9: (Color online) (a) Reduction in OD and AD regions 
of the system 4HD with 77 for uj = 3.0. (b) and (c) show 

the corresponding transition routes from AD and OD states 
to OS state, (d) The emergence of OS regions in the system 
(IB3D with the increase of 77 for N — 200, uj — 3.0 and p — 120 
in ()B3D . 


Even with N = 200 oscillators in system (IB3D . we have 
shown the reduction in AD and OD regions of the system 
in Fig. (9^d). The temporal behavior of the system for 
different values of e and 77 are shown in Figs. HOT a)-(d). 


f*(w<) 
E ij — 


f Xi - coyi - r\xi 
\Vi+ UXi - TiPi 
(1 - 2Sij)Sn 0 

0 — Si2 


H (w j) = wj, 

g( u W) = ( B1 ) 


The eigenvalues corresponding to the trivial equilib¬ 
rium point ( 0 , 0 , 0 , 0 ) of the system are 

M = ^(2-e + r7) 

± yj (e 2 + 77 2 — 4cj 2 ) d= e 2 r] 2 + 4(e 2 — p 2 )u .; 2 . (B2) 

We demonstrate the destabilization of this AD state as 
well as the OD state corresponding to the system m 
with respect to the feedback in Fig. (9^a), where we can 
find that an increase in 77 causes the reduction in AD 
and OD regions of the system. The transition route fol¬ 
lowed by the system, as it transits to oscillatory state, is 
shown in Figs. [9fb) and[9fc). These figures show that the 
oscillations are resurrecting from the AD state via Hopf 
bifurcation as in the previous cases, whereas the resurrec¬ 
tion of oscillations from OD state occurs through saddle 
node bifurcation. 

Recently, AD and OD in the same system with large 
number of oscillators (TV) coupled globally has been seen 
in [58|, whose equation is defined by 

H(w j) = wj , g(u(t)) = -^u(t) 

GX=i <M% + <W)) ) B3) 


_ f i(l -NSa) 
0 


2. Conjugate coupling 

Next, we consider the case of N = 2 oscillators coupled 
through a conjugate coupling described by 

H(w i ) = Wj , L ^ = ( 1 _°^ q ) , g(u(i)) = Iu(t). (B4) 

This system has a trivial equilibrium point at eo: 
(0, 0, 0, 0) and has pairs of non-trivial equilibrium points 

for e > ui, ei, 2 : (aj, bj, -a*, -b*), where aj = 
b\ = —ca\ and e 3 ^: —b^), where a\ = 

zb yj 1+^2 , b% = ca?, in which c = The trivial 

equilibrium point eg has the eigenvalues 


Mi ,2 = (l + 77)±i v /cj(cj + e), 

M3, 4 = 1 ± iy/u(u) e). (B5) 

For 77 = 0, the equilibrium point eo is unstable for all 
values of e and for e = uj a saddle node type bifurcation 
occurs which stabilizes e\ and e 2 as shown in Fig. HH a). 
For 77 7 ^ 0, the eigenvalues corresponding to the equilib¬ 
rium point eo are still unstable, and thus the system is 
still free of AD. Then, the equilibrium points e\ and 
are also destabilized through Hopf bifurcations which is 
demonstrated in Fig. fTTf b). Then, the OD regions of the 
system in the (e, 77 ) space is shown in Fig. QjJc), which 
clearly demonstrates the destabilization of OD with the 
introduction of feedback. 
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Coupling Strength e 


FIG. 11: (Color online) (a) The quenching of oscillations 
through saddle-node bifurcation (SNB) while 77 = 0 in the 
case of conjugate coupling, (b) The destabilization of the OD 
state with the introduction of feedback via Hopf bifurcation 
(. HB ). (c) The reduction in the OD regions of the system in 
the (e, 77 ) space for uj — 5.0. 
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FIG. 13: (Color online) Reduction in OD regions of the sys¬ 
tem (IB71) (a) in the (/i, 77 ) space for e = 1.0, (b) in the (e,rj) 
space for ^ = 1.0. In both the figures, we fixed a = 0.1, 
b — 0 . 1 , c — 18.0 and k — 1 . 0 . 


4. Direct and indirect coupling in Rossler system 

We consider the case of the N = 2 coupled chaotic 
Rossler system [59| defined by 


FIG. 12 : (Color online) (a) shows the AD, OD and OS regions 
of the system (IB6D with respect to e and 77 for b = 0.5. (b) 
The transition route from OD state to OS state via AD state. 


3. Repulsive link: van der Pol oscillator 


Next, we illustrate the role of feedback in the case two 
van der Pol oscillators coupled diffusively through a re¬ 
pulsive link. The corresponding dynamical equations are 
defined through jl 2 | 


f< ( w< )= ( 6(1-4)*-a*)- H (^) = w i- 

= ( t 1 1 iy ) ’ = QU ^’ 

Q=K^)- ( B6 ) 

The AD and OD regions of the system in the (e, 77 ) space 
are given in Fig. [ 12 ] which show the reduction in the AD 
and OD regions with respect to 77 . The transition from 
AD to oscillatory state occurs through a Hopf bifurca¬ 
tion. On the other hand considering the transition from 
OD to oscillatory state, the OD state is transformed to 
AD through inverse pitchfork bifurcation and the oscil¬ 
latory state arises from the AD state through Hopf bi¬ 
furcation. This shows that the role of the feedback in 
setting oscillations back in the system is not restricted to 
any particular oscillator. In the next example, we show 
that the feedback can destabilize oscillation quenching 
scenario even in chaotic oscillators. 


fi(Wi) 

D ij — 


( -y% - zi \ 

Xi + ayi , H(wj)=Wj, 

b + zi{xi - c) ) 


(1-2*,.) + #^ 0 0\ j 

0 0 0, g(u(t)) = -U {t). 

0 0 0 / 



(B7) 


Here, the oscillators are coupled directly by diffusive 
type coupling and are also coupled indirectly to an en¬ 
vironment defined by the variable v. This system (1B71) 
has only non-zero equilibrium points which include e\p. 

(x*,y*,z*,v*) = ( X *,=sl, where, z* = 


( - q,' 1 ± 7 ; . /(// — c ) 2 - 4 Ma+ri) - \ ' p or p ara _ 

metric choice a = 0 . 1 , b = 0 . 1 , c = 18 and k = 1 , the 
equilibrium point e 2 is stable. The OD regions corre¬ 
sponding to the above equilibrium points are given in 
Figs. fl3T a) and fl3T b) with respect to the direct cou¬ 
pling strength (e) and with respect to the indirect cou¬ 
pling strength (/i), respectively. The above figures clearly 
show that the feedback is applicable even for the case of 
chaotic oscillators. 


5. Other chaotic oscillators 

To illustrate further the role of feedback in chaotic 
oscillators, we consider the N = 2 coupled Sprott and 
Lorentz oscillators which are defined respectively by 
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4 6 8 10 

FIG. 14: (Color online) Reduction in the AD and OD regions 

with respect to 77 in the case of (a) the Sprott system (IB 8 D for ^ 

a = 0.225 and (b) the Lorenz system (IB9D for a = 10, 7 = 28 
and b = 2.67. 


(i)Two coupled Sprott systems with repulsive link j57j |. 

/ -ay* \ 

fi(wi) =1 Xi + Zi , H(wj)=w j, 

\Xi+ y\ - Zi J 

/ (1 — 2Si j )S j 2 0 0 \ 

L ij = 0 -Sn 0 , g(u(t)) = -u(t), (B8) 

\ 0 0 0 J z 

and (ii)Diffusive coupling among two Lorenz oscillators 

a, 


FIG. 15: (Color online) The region in which the conditions 
given in 1 IU 211 are satisfied is depicted in the above figure. 
The light-blue shaded region (I) which lies under the light- 
blue curve (curve - 1 ) corresponds to A 3 > 0 , the region under 
the pink curve (curve -2) which is denoted by region (II) cor¬ 
responds to A\ >0, green shaded region (III) corresponds to 
Ao > 0. The regions satisfying the conditions (ii) and (iii) in 
(IC2D are the gray shaded region (IV) (the region under curve 
-4) and red shaded region (V) (the region under curve -5), sep¬ 
arately. The other condition A 2 > 0 is satisfied everywhere 
in the considered region of (e, 77 ). The region in which all 
the conditions in ( 1 C 2 D are satisfied is denoted by blue shaded 
region (VI) (the region under curve - 6 ). 


/ cr{Vi 

~ Xi) 


7 Xi -y%- XiZi 

11 

X 

\ Xiyi 

- bzi 


0 

0 

0 

1 

0 

0 

0 

g(u(*)) = 

(1 - 2 Sij) 

0 

0 

z 


The AD and OD regions corresponding to the coupled 
Sprott and Lorenz systems ( (IB8D and (|B9[0 are given in 
Fig. [TH which again confirm that the feedback wipes out 
AD and OD in the system. 


Aq 2 — J [(4(1 - e) + OJIU 2 )a - (e - 4 )(cj 2 a 2 + u%ai)\ , 
Ao 3 — (1 + ctq)(l + ^ 2 ) — 2 "F UJl ^ 2 )- (01) 

Here a = a\ + < 22 . The eigenvalue equation in (fl3l) can 
be solved directly to get the stable regions of the steady 
state (0,0, 0, 0). On the other hand, we can use the R-H 
criteria to have a closer look at the stable regions of the 
system ©• According the R-H criteria, the stable region 
corresponds to the region in which 


Appendix C: AD region in parametrically 
mismatched system ([2]) 

The coefficients in the characteristic eigenvalue equa¬ 
tion (USD are given by 

A 3 — e — 4 — dr], 

A 2 = ^a 2 i ? 2 + arj(3 - e) + (6 - 3e + cjf + o;|), 

A 1 = An?7 2 + A 12 tj + Ai 3 , 

Ao = Aoi ?7 + A02TJ + A03 5 

where 

An — ^(e-2 )a 2 , 

A 12 -— ( 2 e — 3 )h — ( UJ1CL2 + cu2 a i): 

A 13 — (3e — 4) +-(e — 4 )(cj 2 + cj|), 

Aqi = i [(1 - e)a 2 + (cJia 2 +cj 2 ai) 2 ] , 


(i) Ai> 0 i = 3, 2,1,0 

(ii) A 3 A 2 - Ax > 0 

(iii) A 3 A 2 A! - A\- A\Aq > 0 (C2) 

Now we consider the above criteria one by one and obtain 
the required stable region of the steady state. 

(i)a. A 3 > 0: This condition will be satisfied when 

e — 4 

V < —. (C3) 

This simple condition promises that above the value of 
77 = , the trivial equilibrium point can never be stable 

and thus oscillation can be revoked by an increase in the 
value of 77 . The region in the (e, 77 ) space in which the 
above condition is satisfied is denoted by I in Fig. [15] 
(i)b. A 2 > 0: When e 2 — 3e + 3 — ujf — < 0, if the 

condition (A 2 > 0) is satisfied for 77 = 0 then it will be 
satisfied for all values of 77 . If the condition A 2 > 0 is 
not satisfied while 77 = 0 then the condition will not be 
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satisfied for any value of 77. When e 2 — 3 e+ 3 — oj\ —> 0 , 
the condition A2 > 0 will be satisfied only for the values 
of 77 given by 


77 < 


~Aq2 — \J Aq 2 ~ 4AqiMq3 

2t4 0 i 


(C6) 


2 

77 > - 

e — 3 + J e 2 — 3e + 3 — oj\ — oj\ 

a 

- 

2 

77 < - 

e — 3 — a /e 2 — 3e + 3 — oj\ — cj 2 

a 

* 


(C 4 ) 


This condition A2 > 0 is satisfied in the whole region 
considered in Fig. [ 15 ] 

(i)c. A\ > 0 : When A 2 2 — 4 An A13 < 0, if the condition 
is satisfied for 77 = 0 then it will continue to be satisfied 
for all values of 77, if the condition is not satisfied for 
77 = 0, by the variation of 77 also the condition remains 
to be unsatisfied. When A 2 2 — 4 AnAis > 0 , it will be 
satisfied only for the values of 77 given by 


—A 12 + \J Af 2 — 4 t 4 iit 4 i 3 
2An 

—A 12 — A\ 2 — 4t4ht4i3 

2An 


(C 5 ) 


The region of (e, 77) in which the above condition is satis¬ 
fied is denoted by region (II) which lies under pink curve 
(or curve - 2 ) in Fig. [ 15 ] 

(i)d. Aq > 0 : Similar to the previous conditions, if 
Aq 2 ~ ^AqiAqs > 0, the condition will be satisfied when 


— Aq2 + yj ^ 4 q 2 ~ 4:A 01 Aq3 

2Aqi 


Otherwise it will be satisfied for all values of 77 only if the 
condition is satisfied for 77 = 0 . The region in which the 
above condition is satisfied is shown in Fig. [ 15 ] as region 
III. 

(ii) A3A2 — A\ > 0: We can find that A3A2 — A\ is a 
cubic polynomial in 77, and the region in which the above 
condition is satisfied is given by gray shaded region (IV). 
The real roots of 77 satisfying the equation A3A2 — A\ = 0 
forms the boundary of the region. 

(iii) A 3 A 2 Ai — A\ — A3A0 > 0 : The region satisfy¬ 
ing this condition is shown by red shaded region (V), 
whose boundary is the solution of the quintic equation 
A 3 A 2 A 1 -A 2 - A 2 A 0 = 0 . 

From Fig. [ 15 ] we can find that the region satisfying all 
the criteria given in (|C2[) is the region bounded between 
the curves 


and 


Aq = 0 or 77 = 


— Aq2 lb y/A q 2 — 4AqiAq3 

2 A 0 i 


A 3 A 2 A! - A 2 - A 3 A 0 = 0. (C 7 ) 


The region is denoted by blue shaded region VI in Fig. 
[T 5 l In this region, the equilibrium point (0,0,0,0) is 
found to be stable or AD occurs in the region. 


[11 G. Saxena, A. Prasad, and R. Ramaswamy, Phys. Rep. 
521 205 (2012). 

[2] A. Koseska, E. Volkov, and J. Kurths, Phys. Rep. 531 
173 (2013). 

[3] A. Koseska, E. Volkov, and J. Kurths, Phys. Rev. Lett. 
Ill 024103 (2013). 

[4] T. Banerjee and D. Ghosh, Phys. Rev. E 89 062902 
(2014). 

[5] K. Konishi, Phys. Rev. E 68 067202 (2003). 

[6] W. Zou, D. V. Senthilkumar, J. Duan, and J. Kurths, 
Phys. Rev. E 90 032906 (2014). 

[7] F. M. Atay, Phys. Rev. Lett. 91, 094101 (2003). 

[8] W. Zou, D. V. Senthilkumar, Y. Tang, Y. Wu, J. Lu, and 
J. Kurths, Phys. Rev. E 88 032916 (2013). 

[9] M. Lakshmanan and D.V. Senthilkumar, Dynamics of 
Nonlinear Time-Delay Systems , (Springer, Berlin, 2010). 

[10] A. Prasad, M. Dhamala, B.M. Adhikari, and R. Ra¬ 
maswamy, Phys. Rev. E 81 027201 (2010). 

[11] W. Zou, D. V. Senthilkumar, A. Koseska, and J. Kurths, 
Phys. Rev. E 88 050901 (R) (2013). 

[12] C. R. Hens, P. Pal, S. K. Bhowmick, P. K. Roy, A. Sen 
and S. K. Dana, Phys. Rev. E 89 032901 (2014). 

[13] G. B. Ermentrout and N. Kopell, SIAM J. Appl. Math 
50 125 (1990). 

[14] R. Curtu, Physica D 239 504 (2010). 


[15] M. Dolnik and M. Marek, J. Phys. Chem 92 2452 (1988). 

[16] M. Toiya, V. K. Vanag, and I. R. Epstein, Angew. 
Chem., Int. Ed. 47, 7753 (2008). 

[17] M. Heinrich, T. Dahms, V. Flunkert, S. W. Teitsworth, 
E. Scholl, New J. Phys 12 113030 (2010). 

[18] D. Ruwisch, M. Bode, D. Volkov, and E. Volkov Int. J. 
Bifurcation chaos Appl. Sci. Eng. 9 1969 (1999). 

[19] B. Gallego and P. Cessi, J. Clim. 14 2815 (2001). 

[20] S. Ching, P. L. Purdon, S. Vijayan, N. J. Kopell, and E. 
N. Brown, PNAS, 109 3095 (2012). 

[21] J. Jalife, R. A. Gray, G. E. Morley, and J. M. Davidenko, 
Chaos 8 79 (1998). 

[22] L. Rubchinsky and M. Sushchik, Phys. Rev. E 62 6440 

( 2000 ). 

[23] W. Zou, D. V. Senthilkumar, M. Zhan, and J. Kurths, 
Phys. Rev. Lett 111 014101 (2013). 

[24] S. Kak, Circuits, Systems, and Signal Processing 12 263 
(1993). 

[25] A. Becskei, B. Seraphin, and L. Serrano, Embo. J. 20 
2528 (2001). 

[26] S. Draghici, Int. J. Neural Syst. 8 113 (1997). 

[27] V. K. Chandrasekar, J. H. Sheeba, and M. Lakshmanan, 
Chaos 20 045106 (2010) 

[28] M. G. Rosenblum and A. S. Pikovsky, Phys. Rev. Lett. 
92 114102 (2004). 



















15 


[29] O. V. Popovych, C. Hauptmann, and P. A. Tass, Phys. 
Rev. Lett. 94 164102 (2005). 

[30] G. Yuan, X. Zhang, and Z. Wang, Optik-Int. J. Light 
Electron Opt. 125 1950 (2014). 

[31] G. F. Franklin, J. D. Powell, and A. Emami-Naeini, Feed¬ 
back control of dynamical systems (Addison-Wesley, 3rd 
Ed., UK, 1994). 

[32] J. E. Gough, R. Gohm, and M. Yanagisawa, Phys. Rev. 
A 78 062104 (2008). 

[33] S. Lloyd, Phys. Rev. A 62 022108 (2000). 

[34] S. Little and P. Brown, Ann. N. Y. Acad. Sci 1265 9 

( 2012 ). 

[35] S. Santaniello, G. Fiengo, and L. Glielmo, in Proceedings 
of the IEEE International conference on Control Appli¬ 
cation (CCA 2008) (IEEE, New York, 2008) pp. 666-671. 

[36] J. A. Anderson, An Introduction to Neural Networks , (A 
Bradford Book, 3rd Ed., (1997)). 

[37] Y. Kuromato, Chemical oscillations, Waves and Turbu¬ 
lence , (Springer, Berlin, 1984). 

[38] C. U. Choe, T. Dahms, P. Hovel, and E. Scholl, Phys. 
Rev. E 81 025205(R) (2010). 

[39] T. Aoyagi, Phys. Rev. Lett. 74 4075 (1995). 

[40] S. Uchiyama, Physica A 391 2807 (2012). 

[41] N. Tukhlina and M. Rosenblum, J. Biol. Phys. 34 301 
(2008). 

[42] O. D’Huys, I. Fischer, J. Danckaert, and R. Vicente, 
Phys. Rev. E 83 046223 (2011). 

[43] K. Liidge, Nonlinear Laser Dynamics , (Wiley-VCH, Ger¬ 
many, 2011) 

[44] P. R. Bandyopadhyay and A. M. Helium, Sci. Rep. 4 
6650 (2014). 

[45] E. M. Izhikevich, Dynamical systems in Neuroscience: 


The Geometry of Excitability and Bursting , (MIT Press, 
Cambridge, 2007). 

[46] N. Tukhlina, M. Rosenblum, A. Pikovsky, and J. Kurths, 
Phys. Rev. E 75 011918 (2007). 

[47] O. V. Popovych, C. Hauptmann, and P. A. Tass, Biol. 
Cybern. 95 69 (2006). 

[48] K. Pyragas, O. V. Popovych, and P. A. Tass, Euro. Phys. 
Lett. 80 40002 (2007). 

[49] J. H. Sheeba, V. K. Chandrasekar, and M Lakshmanan, 
Phys. Rev. Lett 103 074101 (2009). 

[50] G. C. Sethia and A. Sen, Phys. Rev. Lett 112 144101 
(2014). 

[51] L. Schmidt and K. Krischer, Phys. Rev. Lett 114 034101 
(2015). 

[52] A. Zakharova, M. Kapeller, and E. Scholl, Phys. Rev. 
Lett. 112 154101 (2014). 

[53] A. Koseska, E. Volkov, and J. Kurths, Europhys. Lett. 
85 28002 (2009). 

[54] A. Koseska, E. Volkov, and J. Kurths, Chaos 20 023132 

( 2010 ). 

[55] M. Lakshmanan, S. Rajasekar, Nonlinear Dynamics: In- 
tegrability, Chaos and Patterns , (Springer-Verlag, Berlin, 
2003) 

[56] K. Bar-Eli, Physica D 14 242 (1985). 

[57] C. R. Hens, O. I. Olusola, P. Pal, and S. K. Dana, Phys. 
Rev. E 88 034902 (2013). 

[58] M. Nandan, C. R. Hens, P. Pal, and S. K. Dana, Chaos 
24 043103 (2014). 

[59] V. Resmi, G. Ambika, and R. E. Amritkar, Phys. Rev. E 
84 046212 (2011). 



